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r-^s ' We review lattice QCD results for glueballs (including a discussion of mixing with 

scalar mesons), hybrid mesons and exotic mesons (such as B^Bs molecules). We 
also discuss string breaking as a mixing between colour flux states and BB states. 



1 Introduction 



The most systematic approach to non-perturbative QCD is via lattice tech- 
niques. Lattice QCD needs as input the quark masses and an overaU scale 
^\ ' (conventionally given by Aqcd)- Then any Green function can be evaluated 

^^ . by taking an average of suitable combinations of the lattice fields in the vacuum 

^^ ' samples. This allows masses to be studied easily and matrix elements (partic- 

ularly those of weak or electromagnetic currents) can be extracted straightfor- 
,S^ ' wardly. Unlike experiment, lattice QCD can vary the quark masses and can 

Mh, also explore different boundary conditions and sources. This allows a wide 

O ' range of studies which can be used to diagnose the health of phenomenological 

QJ , models as well as casting light on experimental data. 

One limitation of the lattice approach to QCD is in exploring hadronic 
decays because the lattice, using Euclidean time, has no concept of asymp- 
totic states. One feasible strategy is to evaluate the mixing between states 
of the same energy - so giving some information on on-shell hadronic decay 
C^ [ amplitudes. 

For comparison with models and for ease of computation, the special case 
of infinitely heavy sea quarks (namely neglect of quark effects in the vacuum: 
the quenched approximation) is often used. We shall also present results from 
including sea quark effects - usually two flavours of degenerate sea quark of 
mass equivalent to strange quarks or heavier. 

The quark model gives a good overall description of hadronic spectra. Here 
I will discuss lattice results for states which go beyond the quark model: glue- 
balls, exotic mesons and hybrid mesons. I also discuss the related phenomenon 
of string breaking. 



2 Glueballs 

Glueballs are defined to be hadronic states made primarily from gluons. The 
full non-pcrturbative gluonic interaction is included in quenched QCD. In the 
quenched approximation, there is no mixing between such glueballs and quark 
- antiquark mesons. A study of the glueball spectrum in quenched QCD is 
thus of great value. This will allow experimental searches to be guided as well 
as providing calibration for models of glueballs. A non-zero glueball mass in 
quenched QCD is the "mass-gap" of QCD. To prove this rigourously is one 
of the major challenges of our times. Here we will explore the situation using 
computational techniques. 

In principle, lattice QCD can study the meson spectrum as the sea quark 
mass is decreased towards experimental values. This will allow the unambigu- 
ous glueball states in the quenched approximation to be tracked as the sea 
quark effects are increased. It may indeed turn out that no meson in the phys- 
ical spectrum is primarily a glueball - all states are mixtures of glue, qq, qqqq, 
etc. We shall later discuss lattice results on the mixing of glueballs and scalar 
mesons (ie qq states). 

In lattice studies, dimensionlcss ratios of quantities are obtained. To ex- 
plore the glueball masses, it is appropriate to combine them with another very 
accurately measured quantity to have a dimensionless observable. Since the 
potential between static quarks is very accurately measured from the lattice, 
it is now conventional to use rg for this comparison. Here ro is implicitly 
defined by r^dV(r)/dr = 1.65 at r = ro where V(r) is the potential energy 
between static quarks which is easy to determine accurately on the lattice. 
Conventionally rg ~ 0.5 fm. 

Theoretical analysis indicates that for Wilson's discretisation of the gauge 
fields in the quenched approximation, the dimensionless ratio mro will differ 
from the continuum limit value by corrections of order a?. Thus in fig. 1 
the mass of the J^^=Q'^^ glueball is plotted versus the lattice spacing a?. 
The straight line then shows the continuum limit obtained by extrapolating to 
a = 0. As can be seen, there is essentially no need for data at even smaller 
a-values to further fix the continuum value. The yaki^shown corresponds to 
TO(0"'"+)ro = 4.33(5). Since several lattice groups El^BEra have measured these 
quantities, it is reassuring to see that the purely lattice observables are in exceL 
lent agreeraent. The publicised difference of quoted to(0"'"+) from UKQCDH 
and GFll eI comes entirely from relating quenched lattice measurements to 
values in GeV. 

In the quenched approximation, different hadronic observables differ from 
experiment by factors of up to 10%. Thus using one quantity or another to 
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set the scale, gives an overall systematic error. Here I choose to set the scale 
by taking the conventional value of the string tension, ^/a = 0.44 GeV, which 
then corresponds to r,^^ — 373 MeV. An overall systematic error of 10% is then 
to be included to any extracted mass. This yields to(0++) = 1611(30)(160) 
MeV where the second error is the systematic scale error. Note that this is 
the glueball mass in the quenched approximation - in the real world significant 
mixing with qq states could modify this value substantially. 

In the Wilson approach, the next lightest glueballs are aB the tensor 
m(2++)ro = 6.0(6) (resulting in to(2++) = 2232(220) (220) MeV) and the pseu- 
doscalar m(2++)ro = 6.0(1.0). Although the Wilson discretisation provides a 
definitive study of the lightest (0"*""*") glueball in the continuum limit, other 
methods are competitive for the determination of the mass of heavier glue- 
balls. Namely, using an improved gauge discretisation which has even smaller 
discretisation errors than the a^ dependence of the Wilson discretisation, so 
allowing a relatively coarse lattice spacing a to be used. To extract mass val- 
ues, one has to explore the time dependence of correlators and for this reason, 
it is optimum to use a relatively small time lattice spacing. Thus an asymmet- 
ric lattice spacing is most appropriate. The results Q are shown in fig. 2 and 



for low lying states are that m{Q++)ro = 4.21(11)(4), m(2++)ro = 5.85(i)(6), 
m(0-+) = 6.33(7)(6) and m(l+-)ro = 7.18(4)(7). Another recent studyS has 
used an improved discretisation based on the perfect action approach (without 
a space-time asymmetry) and obtains results consistent with earlier work. 

One signal of great interest would be a glueball with J^*^ not allowed for 
qq - a spin-exotic gluebaJi|Or oddball - since it would not mix with qq states. 
These states are found BBO to be high lying: considerably above 2to(0"'""'"). 
Thus they are likely to be in a region very difficult to access unambiguously 
by experiment. 

Within the quenched approximation, the glueball states are unmixed with 
qq, qqqq, etc. Furthermore, the qq states have degenerate flavour singlet and 
non-singlet states in the quenched approximation. Once quark loops are al- 
lowed in the vacuum, for the favour-singlet states of any given J^^' , there will 
be mixing between the ss state, the uu + dd state and the glueball. One way 
to explore this is to measure directly the scalar mass eigenstates in a study 
with Nf — 2 flavours of sea-quark. Most studies show no significant change 
of the glueball spectrum as dynamical <|jpf'k effects are adckd - but the sea 
quark masses used are still rather large Era. A recent study eI, however, does 
find evidence for a reduced mass, albeit with a rather large lattice spacing, see 
fig. 1. This effect could be due to mixing of scalar mesons and glueballs, as we 
discuss below, or might just be a sign of an enhanced order a? correction at 
the relatively large lattice spacing used. 

Let us now discus the mixing of the scalar glueball and scalar mesons. The 
mass spectrum of qq states has been determined on a quenched lattice arid 
the scalar mesons are found to lie somewhat lighter than the tensor states El. 
These 2"'""'" mesons are experimentally almost unmixed and so will be quite 
close to the quenched mass determination. This suggests that the quenched 
scalar masses from the lattice are at around_L2. GeV and 1.5 GeV (for nfi 
and ss respectively). An independent study HO suggests that the scalar ss 
state is about 200 MeV lighter than the glueball which is a broadly compatible 
conclusion. Thus the glueball, at around 1.6 GeV, lies heavier than the lightest 
qq scalar states. This information can then be combined with mixing strengths 
to give the resulting scalar spectrum. 

It is possible to measure the mixing strength on a quenched lattice even 
though no mixing actually occurS|_Qn a rather coarse lattice (a"-'^ w 1.2 GeV), 
two groups have attempted thislid'u. Their results expressed as the mixing 
for two degenerate quarks of mass arouiid the strange quark mass are similar, 
namely E « 0.3 GeVEiJ and 0.5 GeVB. From this evaluation of the mixing 
strength, one can use a mass matrix to estimate the mass shift induced in the 
glueball and scalar meson. The relevant mass matrix is (in a glueball, qq basis 
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in GeV units): 
1.2 0.5' 
0.5 1.8^ 

which would give a downward shift of the glueball mass by 20%. This is in 
qualitative agreement with our direct determination with Nj = 2 flavours of 
sea-quark that the lightest scalar mass is reduced significantly at this lattice 
spacing as shown in fig. 1. 

Note that at this coarse lattice spacing the quenched glueball mass is re- 
duced (see fig. 1) below the canonical value of 1.6 GeV. Thus a study at smaller 
lattice, spacing is needed. An exploratory attempt to extrapolate to the contirtj 
uumliil gave a very small mixing of 86(64) MeV, while the other determinationQ 
uses clover improvement so order a effects in the extrapolation to the contin- 
uum are suppressed and one would not expect a significant decrease in going 
to the continuum limit. What this discussion shows is that precision studies 
of the mixing on a lattice have not yet been achieved. 

As well as this mixing of the glueball with qq states, there will be mixing 
with qqqq states whkJi will be responsible for the hadronic decays. A first 
attempt to study thisllj yields an estimated width for decay to two pseudoscalar 
mesons from the scalar glueball of order 100 MeV. A more realistic study would 
involve taking account of mixing with the nn and ss scalar mesons as well. 



3 Exotic states 

By exotic state we mean any state which is not dominantly a qq or qqq state. 
Such examples have been known for a long time: the deuteron is a proton- 
neutron molecule for example. It is very weakly bound (2 MeV) and is quite 
extended. Similar molecular states involving two mesons have been conjec- 
tured. 

One case which is relatively easy to study is the BB system, idealised as 
two static quarks and two light quarks. Then a potential as a function of the 
separation R between the static quarks can be determined. Because the static 
quark spin is irrelevant, th£_^states can be classified by the light quark spin 
and isospin. Lattice results E^ (using a light quark mass close to strange) have 
been obtained for the potential energy for Iq — 0,1 and Sq = 0, 1. For very 
heavy quarks, a potential below 2Mb will imply binding of the BB molecules 
with these quantum numbers and L — 0. For the physically relevant case of 
b quarks of around 5 GcV, the kinetic energy will not be negligible and the 
binding energy of the BB molecular states is less clear cut. One way to estimate 
the kinetic energy for the BB case with reduced mass circa 2.5 GeV is to use 
analytic approximations to the potentials found. For example the /g, Sq={0,0) 
case shows a deep binding aX R = which can be approximated as a Coulomb 
potential of — 0.1/i? in GeV units. This will give a di-meson binding energy 
of only 10 MeV. For the other interesting case, {Iq, Sq)~{0,l), a harmonic 
oscillator potential in the radial coordinate of form — 0.04[1— (r — 3)^/4] in GeV 
units leads to a kinetic energy which completely cancels the potential energy 
minimum, leaving zero binding. This harmonic oscillator approximation lies 
above the estimate of the potential, so again we expect weak binding of the 
di-meson system. 

Because of these very small values for the di-meson binding energies, one 
needs to retain corrections to the heavy quark approximation to make more 
definite predictions, since these corrections are known to be of magnitude 46 
MeV from the B, B* splitting. It will also be necessary to extrapolate the 
light quark mass from strange to the lighter w, d values to make more definite 
predictions about the binding of BB molecules. 

Models for the binding of two B mesons_involve, as in the case of the 
deuteron, pion exchange. The lattice study Ej is able to make a quantita- 
tive comparison of lattice pion exchange with the data described above and 
excellent agreement is obtained at larger R values, as expected. 



4 Hybrid Mesons 

By hybrid meson, I mean a meson in wliicli tlie gluonic degrees of freedom are 
excited non-trivially. I first discuss hybrid mesons with static heavy quarks 
where the description can be thought of as an excited colour string. I then 
summarise the situation concerning hght quark hybrid mesons. 

Consider QQ states with static quarks in which the gluonic contribution 
may be excited. We classify the gluonic fields according to the symmetries of 
the system. This discussion is very similar to the description of electron wave 
functions in diatomic molecules. The symmetries are (i) rotation around the 
separation axis z with representations labelled by Jz (ii) CP with represen- 
tations labelled by ^(-l-) and u(— ) and (iii) C7?.. Here C interchanges Q and 
Q, P is parity and 7?, is a rotation of 180*^ about the mid-point around the y 
axis. The C7?. operation is only relevant to classify states with J^ = 0. The 
convention is to label states of J^ = 0, 1, 2 by S, II, A respectively. The ground 
state (E+ — g) will have Jz — ^ and CP — -f-. 

The exploration of tjifi_epergy levels of other representations has a long 
history in lattice studies tSlla. The first excited state is found to be the n„. 
This can be visualised as the symmetry of a string bowed out in the x direction 
minus the same deflection in the —x direction (plus another component of the 
two-dimensional representation with the transverse direction x replaced by y) , 
corresponding to flux states from a lattice operator which is the difference of 
U-shaped paths from quark to antiquark of the form Fl — U. 

Recent lattice studies lj have used an asymmetric space/time spacing 
which enables excited states to be determined comprehensively. These re- 
sults confirm the finding that the n„ excitation is the lowest lying and hence 
of most relevance to spectroscopy. 

From the potential corresponding to these excited gluonic states, one can 
determine the spectrum of hybrid quarkonia using the Schrodinger equation 
in the Born-Oppenheimer approximation. This approximation will be good if 
the heavy quarks move very little in the time it takes for the potential between 
them to become established. More quantitatively, we require that the potential 
energy of gluonic excitation is much largcrjthan the typical energy of orbital or 
radial excitation. This is indeed the caseEJ, especially for 6 quarks. Another 
nice feature of this approach is that the self energy of the static sources cancels 
in the energy difference between this hybrid state and the QQ states. Thus the 
lattice approach gives directly the excitation energy of each gluonic excitation. 

The IIu symmetry state corresponds to excitations of the gluonic field in 
quarkonium called magnetic (with L^'-^ = 1"' ) and pseudo-electric (with 1 ^) 
in contrast to the usual P-wave orbital excitation which has L^'~^ = 1^^. Thus 



we expect different quantum number assignments from those of the gluonic 
ground state. Indeed combining with the heavy quark spins, we get a degener- 
ate set of 8 states with J^^ = 1~, 0"+, 1"+, 2"+ and 1++, 0+", 1+", 2+" 
respectively. Note that of these, J^*^ = 1 '", 0^ and 2^ are spin-exotic 
and hence will not mix with QQ states. They thus form a very attractive goal 
for experimental searches for hybrid mesons. 

The eightfold degeneracy of the static approach will be broken by various 
corrections. As an example, one of the eight degenerate hybrid states is a 
pseudoscalar with the heavy quarks in a spin triplet. This has the same overall 
quantum numbers as the S-wave QQ state {rjb) which, however, has the heavy 
quarks in a spin singlet. So any mixing between these states must be mediated 
by spin dependent interactions. These spin dependent interactions will be 
smaller for heavier quarks. It is of interest to establish the strength of these 
effects for b and c quarks. Another topic of interest is the splitting between the 
spin exotic hybrids which will come from the different energies of the magnetic 
and pseudo-electric gluonic excitations. 

One way to go beyond the static approach is to use the NRQCD approx- 
imatioii. which then enables the spin dependent effects to be explored. One 
studyll3 finds that the L^'-^ = 1^ and 1 ^ excitations have no statistically 
significant splitting although the 1^ excitation does lie a little lighter. This 
would imply, after adding in heavy quark spin, that the J^*^ = 1 '" hybrid 
was the lightest spin exotic. Also a relatively large spin splitting was foundliZl 
among the triplet states considering, however, only magnetic gluonic excita- 
tions. 

Confirmation of the ordering of the_smn exotic states also comes from 
lattice studies with propagating quarkslj'ESIlj which are able to measure masses 
for all 8 states. We discuss this evidence in more detail below. 

Within the quenched approximation, the lattice evidence for bb quarks 
points to a lightest hybrid spin exotic with J^'-' = 1 ^ at aii energy gpen by 
{■mH~m2s)j-a =1.8 (static potential^); 1.9 (static potentiafi, NRQCD0); 2.0 
(NRQCDEj). These results can be summarised as {ruH — rn2s)i'o = 1.9 ± 0.1. 
Using the experimental mass of the T(2S'), this implies that the lightest spin 
exotic hybrid is at rriH = 10.73(7) GeV including a 10% scale error. Above 
this energy there will be many more hybrid states, many of wiiich will be spin 
exotic. A discussion of hybrid decay channels has been given til. 

The excited gluonic static potential has also been determined including 
sea quarks (Nj = 2 flavours) and no significant difference is seenE3- see fig. |[ 
Thus the quenched estimates given above are not superseded. 

I now focus on lattice results for hybrid mesons made from light quarks 
using fully relativistic propagating quarks. There will be no mixing with qq 
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Figure 3: The potential energjp^pr quenched and 2 flavours of sea quark for the ground state 
and first excited gluonic stateEj- 



mesons for spin-exotic hybrid mesons and these are of special interest. The 
first study of this area was by the UKQCD CollaborationEl who used operators 
motivated by the heavy quark studies referred to above to study all 8 J^'-^ 
values coming from iJ' = 1"' and 1 '" excitations. The resulting mass 
spectrum gives the J^*^ = 1 ^ state as the lightest spin-exotic state. Taking 
account of the systematic scale errors in the lattice determination, a mass of 
2000(200) MeV is quoted for this hybrid meson with ss light quarks. Although 
not directly measured, the corresponding light quark hybrid meson would be 
expected to be around 120 MeV lighter. 

A second lattice group has also evaluated hybrid meson spectra with prop- 
agating quarks from quenched lattices. They obtain E3 masses of the 1 '" 
state with statistical and various systematic errors of 1970(90) (300) MeV, 
2170(80) (100) (100) MeV and 4390(80) (200) MeV for nn, ss and cc quarks 
respectively. For the 0^' spin-exotic state they have a noisier signal but ev- 
idence that it is heavier. They also explore mixing matrix elements between 
spin-exotic hybrid states and 4 quark operators. 

A first attempt has been madeEj to determine the hybrid meson spectrum 
using full QCD. The sea quarks used have several different masses and an 
extrapolation is made to the limit of physical sea quark masses, yielding a 
mass of 1.9(2) GeV for the lightest spin-exotic hybrid meson, which again is 



found to be the 1 ^. In principle this calculation should take account of sea 
quark effects such as the mixing between such a hybrid meson and qqqq states 
such as rjTT, although it is possible that the sea quark masses used are not light 
enough to explore these features. 

The three independent lattice calculations of the light hybrid spectrum are 
in good agreement with each other. They imply that the natural energy range 
for spin-exotic hybrid mesons is around 1.9 GeV. The J^*^ = 1 '" state is found 
to be lightest. It is not easy to reconcile these lattice results with experimental 
indications Eil for resonances at 1.4 GeV and 1.6 GeV, especially the lower 
mass value. Mixing with qqqq states such as rjn is not included for realistic 
quark masses in the lattice calculations. This can be interpreted, dependent on 
one's viewpoint, as either that the lattice calculations are incomplete or as an 
indication that the experimental states may have an important meson-meson 
component in them. 

5 String breaking 

Although not directly related to exotic states as such, a discussion of string 
breaking is relevant to a study of the confinement mechanism. Here I concen- 
trate on the situation of a static quark and antiquark at separation R with a 
gluonic field in a given representation labelled by Jz and CP where the separa- 
tion axis is z. The ground state gluonic configuration has Jz = and CP = + 
(also CTZ — +, see discussion above). 

If a light quark and antiquark are to be produced, they must respect these 
symmetries. When these light quarks have no angular momentum about the 
separation axis, they will either be in a singlet state {Sz — 0, CP = — ) or a 
triplet state {Sz = or S'^ = 1, CP = +). Thus the only possible transition is 
to the 5*2 = triplet state. This is known as the ^Pq model, although here it is 
exact. Having established the allowed transitions, I now discuss their strength. 

Consider the energy of a static QQ versus R for two different situations: 
(i) a gluonic flux in the ground state (this is the static potential V{R)) and 
(ii) the energy of two ground state mesons {Qq plus Qq) at separation R (ie 
as B and B mesons but with static heavy quark). We expect (i) to be lowest 
lying at small R and (ii) to be lowest lying at large R. Where these energies 
become degenerate is the string-breaking separation Rb- This is around 
fm - see fig. y. As has been emphasised in a study of adjoint string breakingE 
the lattice enables string breaking to be studied as a mixing phenomenon. 

In the quenched approximation there will be no string breaking but the 
matrix element can be estimated. In full QCD studies with sea quarks, string 
breaking is enabled and there will be mixing between these states. If this 
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mixing were substantial, the "potential energy" determined from Wilson loops 
should show a saturation at energy V{R) = 2mQq for R > Kb- Unfortunately, 
this effect has not yet been seen directly in lattice studies. 

A comprehensive lattice study l3 has been made with static quarks and 
light quarks of mass around the strange mass. This study finds that observa- 
tion of the direct splitting at i? « Rb has errors too large to be conclusive. 
The magnitude of the mixing matrix element can, however, be evaluated with 
adequate accuracy. This suggests that indeed the mixing is weak (energy mix- 
ing strength of a few tens of MeV only). Such a weak mixing would explain 
why no direct evidence has been seen for string breaking from lattice studies 
of Wilson loops. 

This discussion of string breaking can be extendedEJ to cover hybrid meson 
decays as well. 

6 Conclusions 

Quenched lattice QCD is well understood and accurate predictions in the con- 
tinuum limit are increasingly becoming available. The lightest glueball is scalar 
with mass m(0++) — 1611(30)(160) MeV where the second error is an over- 
all scale error. The excited glueball spectrum is known too. The quenched 
approximation also gives information on quark-antiquark scalar mesons and 
their mixing with glueballs. This determination of the mixing in the quenched 
approximation also sheds light on results for the spectrum directly in full QCD 
where the mixing will be enabled. There is also some lattice information on 
the hadronic decay amplitudes of glueballs. 

Evidence exists for a possible BsBs molecular state. 

For hybrid mesons, there will be no mixing with qq for spin-exotic states 
and these are the most useful predictions. The J^^ = 1 '" state is expected 
at 10.73(7) GeV for h quarks, 2.0(2) GeV for s quarks and 1.9(2) GeV for 
M, d quarks. Mixing of spin-exotic hybrids with qqqq or equivalently with 
meson-meson is allowed and will modify the predictions from the quenched 
approximation. 

String breaking at a separation of 1.2 fm has been studied and the breaking 
matrix element is found to be weak. 
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